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Abstract 

The mathematical apparatus of quantum-mechanical angular mo- 
mentum (re)coupling, developed originally to describe spectroscopic 
phenomena in atomic, molecular, optical and nuclear physics, is em- 
bedded in modern algebraic settings which emphasize the underlying 
combinatorial aspects. SU(2) recoupling theory, involving Wigner's 
3nj symbols, as well as the related problems of their calculations, 
general properties, asymptotic limits for large entries, plays nowadays 
a prominent role also in quantum gravity and quantum computing 
applications. We refer to the ingredients of this theory -and of its 
extension to other Lie and quantum group- by using the collective 
term of 'spin networks'. 

Recent progress is recorded about the already established connec- 
tions with the mathematical theory of discrete orthogonal polynomi- 
als (the so-called Askey Scheme), providing powerful tools based on 
asymptotic expansions, which correspond on the physical side to var- 
ious levels of semi-classical limits. These results are useful not only 
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in theoretical molecular physics but also in motivating algorithms for 
the computationally demanding problems of molecular dynamics and 
chemical reaction theory, where large angular momenta are typically 
involved. As for quantum chemistry applications of these techniques 
include selection and classification of complete orthogonal basis sets 
in atomic and molecular problems, either in configuration space (Stur- 
mian orbitals) or in momentum space. In this paper we list and discuss 
some aspects of these developments -such as for instance the hyper- 
quantization algorithm- as well as a few applications to quantum grav- 
ity and topology, thus providing evidence of a unifying background 
structure. 

1 Introduction 

The (re)coupling theory of many SU{2) angular momenta -the Racah-Wigner ten- 
sor algebra- is the most exhaustive framework in dealing with interacting many- 
angular momenta quantum systems that can be modeled by means of conservative 
polylocal two-body interactions pp. The essential features of this algebra can be 
encoded, for each fixed number N = (ra+ 1) of angular momentum variables, into a 
combinatorial object, the spin network graph. Vertices are associated with finite- 
dimensional, binary coupled Hilbert spaces while edges correspond to either phase 
or Racah transforms (implemented by 6j symbols) acting on states in such a way 
that the quantum transition amplitude between any pair of vertices is provided by 
a suitable 3nj symbol. 

One of the goals of this paper is to promote this combinatorial setting to 
(families of) computational quantum graphs (SecGJ), to be thought of as a sort of 
'abacus' encoding diagrammatical rules encountered in quantum collision theory 
[2] or, more in general, as computational spaces of simulators able to support quan- 
tum algorithms for computational problems arising in mathematics and theoretical 
physics [3]. Progress in the semiclassical limits is accounted for in Sec. A few 
examples are then briefly addressed, emphasis being placed on those semiclassical 
versions of the construction, which turn out to play a crucial role in many different 
contexts, ranging from molecular physics (Sec. [4]) and quantum computing and 
discrete quantum gravity models (Sec. [5]). An outlook to applications and a list 
of perspectives for future work are collected in the concluding Sec. El 
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2 $77(2) recoupling schemes as spin network 
graphs 

Important ingredients of modern theoretical physics and chemistry are intimately 
linked to the mathematical theory of classical (as opposed to g-deformed) orthog- 
onal polynomials, in particular of hyper geometric families. The recently proposed 
formalization of the latter within the Askey and Nikiforov schemes [U [5] places 
Racah polynomials -and their g-deformed analogs [6]- at the top of a hierarchy 
from which all the most relevant families are obtained as particular cases or via 
suitable limiting procedures. Although the study of such morphology may be 
considered as an exercise in special function theory, relationships among families, 
addition formulas, linearization formulas and sum rules look like obscure manipu- 
lations of abstract quantities unless a coherent interpretation arising from physical 
applications could be disclosed. Physical applications, on the other hand, can lead 
to new insights into previously unnoticed properties and relationships. 

Dealing with the quantum theory of angular momentum, Racah found a finite 
sum expression for the basic 'recoupling' coefficients or 6j symbols, representing 
matrix elements for the orthogonal basis transformation between two alternative 
binary coupling schemes of three angular momenta. Extensive developments of 
graphical techniques, initiated by the Yutsis school [7], allow one to deal efficiently 
with recouplings of (n + 1) angular momenta, involving the so-called 3nj symbols 
(see the handbook [8] and the review [T] (Topic 12) also for a complete list of 
references) . 

An introduction and some nomenclature are exhibited in [2], where the 'spin 
networks', as well as the underlying 'moves', are shown to be closely related to the 
tree-like graphical representation of hyperspherical coordinates and harmonics, as 
originally suggested in Ref. [9]. 

Explicit formulas and basic relationships among recoupling coefficients and 
harmonic superpositions (or transplant coefficients) and Racah polynomials can 
be found in [10J and in references therein. The close relationship of both vector 
recoupling coefficients and superposition matrix elements between alternative hy- 
perspherical harmonics with orthogonal polynomials of a discrete variable made 
it possible to develop a common classification scheme, from which a number of 
relations and properties can be derived [2]. The important case of S 3 harmonics, 
discussed in [TTJ from the viewpoint of hydrogenoid orbitals in momentum space, 
is reviewed in [12], where further applications for quantum chemistry are outlined. 

The underlying mathematical background of the developments discussed in [2] , 
involves objects like the 3nj symbols of angular momentum theory (and their next 
of kin, the superposition coefficients) which are candidates for further study as 
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orthogonal polynomials of (n — 1) discrete variables. 

The basic features of spin network graphs are illustrated in Fig. [1] for the 
significant case of binary coupling schemes of four angular momenta and their 
recoupling represented by a Wigner 9j coefficient. We refer to [2] for more details 
on their use as 'abacus', namely as graphical devices for practical calculations. For 
the construction of the abacus a proper look in terms of 'presentation' of discrete 
groups shows that the basic tools are the three operations corresponding to the 
presentation of the icosahedral group [131 114j . From a geometrical viewpoint, the 
abacus essentially mixes (according to a well-known presentation of the icosahedral 
symmetries) pentagonal and hexagonal cycles, whose deep origin can be traced in 
abstract algebraic theories [15] . 



According to Bohr correspondence principle, classical concepts become increasingly 
valid in regimes where all (or just a few) quantum numbers are 'large' (as will be 
illustrated in Sec. such regimes are quite commonly encountered in every-day 
analysis of atomic and molecular dynamical processes). 

In handling with angular momenta variables measured in units of h, the clas- 
sical limit H — ► implies in particular that, for finite angular momenta, the j— 
quantum numbers are much larger than one. 

In the case in which all the six angular momenta are large, the Racah 6j sym- 
bol admits a well defined asymptotic limit, whose absolute square (probability) 
corresponds to the classical limit of the associated physical observables. Following 
[161 El H] an d p] (Topic 9) (where a self contained treatment of various asymp- 
totics of angular momentum functions is given, together with the list of original 
references), recall that the square of the symbol has the limiting value given by 
the Wigner formula 



where V is the Euclidean volume of the tetrahedron formed by the six angular 
momentum 'vectors'. 

A major breakthrough in semiclassical analysis is provided by the Ponzano- 
Regge asymptotic formula for the 6j symbol [16J 



where the limit is taken for all entries 3> 1 (recall that K = 1) and t r = j r + 1/2 
with {j r } = {a,b,c,d,e,f}. V is now the Euclidean volume of the tetrahedron 



3 Semiclassical spin networks 
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Figure 1: Spin networks are illustrated here for the case of two quantum systems, 1 and 
2, each endowed with both spin and orbital angular momentum, giving a state with total 
angular momentum j (coupling described by spin-orbit terms in the Hamiltonian) . The 
two alternative ways of pairing s\,S2,li,li, traditionally referred to as Is and jj couplings, 
are encoded in the top left and right 'binary trees'. Their recoupling -achieved by joining- 
edges with identical labels- leads to the 3-valent, not planar graph (of the Yutsis-type, 
as referred to in the literature). The associated recoupling matrix elements give the 9j 
symbol, represented by the 3x3 array on the left-hand side of the formula. In the 
center a portion of what we call an 'abacus', namely an aid in obtaining the procedure 
for getting the decomposition illustrated in the rest of the external region of the figure: 
the three complete graphs on 4 vertices are associated with suitably labelled 6j symbols, 
while dotted arrows correspond to phase changes. The structure of the spin network graph 
displays continuous and dashed lines as well as typical pentagonal and hexagonal cycles 
(providing an overall 'fullerene'-like pattern), associated with the Biedenharn-Elliott and 
Racah identities, respectively. 
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with edges of lengths {£ r }, calculated by using the Cayley determinant (note the 
shift j — > j + 1/2 with respect to the variables employed in calculating the volume 
in (pQ)) and finally 6 r is the angle between the outer normals to the faces which 
share the edge i r . 

The probability amplitude (|2|) has the form of a semiclassical wave function 
since the factor \j\f 7 2AnV is slowly varying with respect to the spin variables, while 
the exponential is a rapidly oscillating dynamical phase (such behavior complies 
with the fact that the square of this asymptotic expression reproduces Wigner's 
result (pQ)). Moreover, according to Feynman path sum interpretation of quantum 
mechanics, the argument of the exponential represents a classical action, and in- 
deed it can be read as P Q f° r pairs (p, q) of canonical variables (angular momenta 
and conjugate angles). A mathematically sophisticated analysis of Ponzano-Regge 
result can be found in [19], where recent developments in the study of the asymp- 
totics of 3j and 6j symbols are also reviewed. 

An interesting issue -arising in connection with the interpretation of spin net- 
works as computational quantum graphs- concerns the phenomenon of disentagle- 
ment. One of the main features of states belonging to the binary coupled Hilbert 
spaces, whose 'recoupling' gives the 6j Racah coefficient, is to represent effectively 
'entangled' quantum states (namely states that cannot be reduced to a product 
of states containing quantum numbers of the individual components). Then the 
Racah transform (implemented by means of a 6j symbol, see the lefthand side of 
([2])), being a quantum transition amplitude, takes care of the fact that the oper- 
ators and Jg, whose eigenvalues are d(d + 1) and e(e + 1) respectively, do not 
belong to the same set of mutually commuting operators and thus cannot be mea- 
sured simultaneously. On the other hand, when reaching the semiclassical limit as 
can be seen from the righthand side of ([2]), the six entries of the 6j symbol appear 
to be on the same footing, a feature that can be viewed as a 'disentanglement' of 
the underlying 'semiclassical' spin networks. Work is in progess for what concerns 
the 9j symbol [20]. 

4 Molecular Physics 

4.1 Atomic collisions and molecular spectroscopy 

A presentation of angular momentum theory from the viewpoint of these appli- 
cations, is given by Zare [21]. The theory developed about thirty years ago in 
|22j dealt with five alternative representations for the quantum-mechanical close- 
coupling formulation [22, 23, 24J of the motion along the internuclear distance of a 
vibrating diatomic molecule or colliding atoms having internal (spin and/or elec- 
tronic) angular momenta. This unified frame transformation approaches of atomic 
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collision theory and concepts of molecular spectroscopy, is originally due to Hund. 
The physical picture and the relevant nomenclature are reviewed in Ref [25J, see 
also [26]. Explicitly |27j . starting from a sum rule equivalent to the Biedenharn- 
Elliott relation, (a pentagonal closed path on the abacus of Fig. Q]) we obtain the 
definition of a 6j symbol as a sum of four 3j symbols by taking a proper limit, 
since one angular momentum is much larger than the others [25J. 

The coupling schemes of four angular momenta are illustrated in Ref. [25] as 
the basic ingredients underlying the classification of the five Hund cases and the 
relationships among representations. It was shown that the transition from one 
coupling scheme to another is performed by an orthogonal transformation whose 
matrix elements can be written in terms of 6j symbols. In the pentagonal arrange- 
ment of five alternative coupling schemes for four angular momenta (represented 
by the tree-like graphs at the vertices) [27], connections (the sides of the pentagon) 
are realized by orthogonal matrices and are related to 6j symbols. This was the 
archetype for the abacus in Fig. (TJ Recent extensions of this approach have shown 
that in a general theory of interacting open shell atoms, 3nj symbols up to n = 6 
occur [28] . 

4.2 Reaction dynamics: Hyperquantization 

An important message that we have learned from work reviewed in Sec. 14. H is the 
view that a continuous variable limit is obtained at high angular momenta from 
the discrete structure typical of quantum mechanics. The opposite viewpoint can 
be considered as well, namely the semiclassical limit may describe a continuous 
structure, and quantum angular momentum algebra provides discretization. 

The search for both alternative reference frames and angular momentum cou- 
pling schemes has been a major challenge in quantum mechanics since its origin, 
and tranformations among them are represented by vector-coupling and recoupling 
coefficients, respectively. The relevant equations can then be formulated in terms 
of quantum numbers, which approximately correspond to constants of motion of 
the systems under study. Fundamental advances have been achieved over the years 
within this framework: in the last Fifties Jacob and Wick introduced the helic- 
ity formalism, widely used for the theoretical treatment of a variety of collisional 
problems; extending Hund's introduction of alternative coupling schemes for a di- 
atomic molecule carrying electronic, spin, and rotational angular momenta (See 
Sec. 14.1ft . These developments fit into the frame transformation theory pioneered 
by Fano and coworkers in the Seventies. 

Indeed, for general anisotropic interaction, discretization procedures can be 
introduced by exploiting Racah algebra, which fosters the introduction of alter- 
native coupling schemes labeled by 'artificial' quantum numbers. This method 
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has been shown to provide an elegant and powerful tool for the solution of the 
reactive scattering Schrondinger equation and, at present, a considerable number 
of methods have been developed in this spirit, among which the 'hyperquantiza- 
tion' algorithm, outlined in greater detail in a number of references |29t [271 130J . 
The technique relies on the hyperspherical coordinate approach when used for 
few-body processes, including rearrangement. For instance, in a reactive triatomic 
process, the reaction coordinate is represented by the hyperradius, which is a mea- 
sure of the total inertia of the system, and an adiabatic representation of the total 
eigenfunction with respect to this coordinate is adopted. In such a way, a quan- 
tization problem on the surface of the a hypersphere (in this case the sphere in a 
6-dimensional Euclidean space) must be solved for fixed values of the hyperradius. 
Then coupled-channel equations are integrated applying a standard propagation 
procedure. The success of this approach is strictly dependent on the accuracy 
and the effectiveness of the method used to solve the fixed hyperradius problem. 
The computation of the adiabatic eigenvalues containing detailed information on 
the structure, rotations, and internal modes parametrically in the hyperrradius is 
typically very demanding. The hyperquantization algorithm exploits the peculiar 
properties of the discrete analogues of hyperspherical harmonics, i.e. generalized 
3j symbols or Hahn polynomials, orthogonal on a grid of points that span the in- 
teraction region [31]. The computationally advantageous aspect of this algorithm, 
besides the elegance of unifying under the language of angular momentum theory 
the dynamical treatment of a reaction, is the structure of the Hamiltonian matrix: 
its kinetic part is simple, universal, highly symmetric, and sparse, while the poten- 
tial displays the diagonal form characteristic of the stereodirected representations 
of the previous section. The hyperquantization algorithm, when implemented for 
reactive scattering calculations |27j , allows considerable savings in memory require- 
ments for storage and in computing time for the building up and diagonalization 
of large basis sets, exploiting the sparseness and the symmetry properties of the 
Hamiltonian matrix. 



5 Quantum gravity and quantum computing 

5.1 Discretized quantum gravity in dimension 3 

There exists an intriguing physical interpretation of the Ponzano-Regge asymp- 
totics ([2]) once we recognize that the expression in the exponential represents the 
classical Regge action |32j -namely the discretized version of Einstein-Hilbert ac- 
tion of General Relativity - for the tetrahedron associated with the 6j symbol in 
the semiclassical limit. 

In Fig. Q] the tetrahedral symmetry of each 6j symbol is encoded into the 
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corresponding 4-vertex graph, where the six entries of the symbol are associated 
with edges and its four triads with faces of a tetrahedron embedded in Euclidean 
3-space (see e.g. [8] for more details on the symmetries and the geometrical 
properties of this symbol). 

In Regge's approach to General Relativity, the edge lengths of a 'triangulated 
spacetime' are taken as discrete counterparts of the metric tensor appearing in 
the usual action for gravity and angular variables (deficit angles) are related to 
the scalar curvature obtained from the Riemann tensor. Technically speaking, a 
Regge spacetime is a piecewise linear (PL) manifold of dimension D dissected into 
simplices, namely triangles in D = 2, tetrahedra in D = 3, 4-simplices in D = 4 
and so on. 'Regge Calculus' became in the early 80's the starting point for a 
novel approach to quantization of General Relativity known as simplicial quantum 
gravity (see [33J and references therein). The quantization procedure most com- 
monly adopted is the Euclidean path sum approach, namely the discretized version 
Feynman path integral describing D-dimensional geometries undergoing 'quantum 
fluctuations'. According to this prescription, the asymptotic functional ([2]) -to be 
understood here as the semiclassical limit of a sum over truly 'quantum' fluctua- 
tions - turns out to be associated with the simplest 3-dimensional 'spacetime', the 
Euclidean tetrahedron. The construction of the so-called Ponzano-Regge 'state 
sum' representing the quantum partition function of simplicial 3-gravity will not 
be discussed further (see e.g. [3; Sec. 5, for an account). 

Several years after Ponzano-Regge paper, a regularized version of their state 
sum -based on representation theory of a 'quantum' deformation of the group 
SU(2)~ was proposed in [33] and shown to be a well-defined (finite) topological 
invariant for closed 3-manifolds. Since then there has been a renewed interest also 
in the asymptotics ([2]) both in connection with the study of 3-manifold geometry 
(and higher-dimensional generalizations [35] ) and in addressing 'loop quantum 
gravity' models, see e.g. [36] and references therein. 




5.2 Quantum automata and topological invariants 

In the past few years there has been a tumultuous activity aimed at introducing 
novel conceptual schemes for quantum computing. The model of quantum simu- 
lator proposed in [371 E] and further discussed in [381 EH] relies on the recoupling 
theory of SU(2) angular momenta and can be viewed as a generalization to ar- 
bitrary values of the spin variables of the usual quantum-circuit model based on 
'qubits' and Boolean gates [3D]. The basic ingredient of such general scheme for 
universal quantum computing are indeed encoded into spin network computational 
graphs of the type depicted in the central portion of Fig. [JJ). Such pictorial rep- 
resentation makes it clear that the computational space of the simulator complies 
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with the architecture of (families of) 'automata^]- According with this kind of in- 
terpretation, a computational process on the spin network can be associated with 
a directed path, namely an ordered sequence of vertices and edges starting from an 
initial quantum state, say \s >j n , and ending in some set of final states {\s > fin}- 
In a series of papers [3H H2J H3J HU H5] families of automata arising from 
the g-deformed analog of the spin network simulator have been implemented in 
order to deal with classes of computationally-hard problems in geometric topol- 
ogy (topological invariants associated with knots and with closed 3-dimensional 
manifolds jj. 

From the point of view of classical complexity theory, computing such invari- 
ants is 'hard', namely could be achieved by a classical computer only by resorting 
to an exponential amount of resources. A computational process which requires an 
amount of resources that grows at most polynomially with the size of the compu- 
tational problem is referred to as 'efficient' (c/r. [45J and references therein for an 
account of algorithmic questions involving braid group and topological invariants 
ok knots). 

In [41[ I42j. efficient (i.e. running in polynomial time) quantum algorithms for 
approximating, within an arbitrarily small range, SU (2) q -colored polynomial in- 
variants of knots have been explicitly worked out. In [33] such algorithms have been 
generalized to deal with 3-manifold invariants, while in [33] connections among 
quantized geometries, topological quantum field theory and quantum computing 
are discussed in detail. 

6 Extensions and applications 

Besides ^-deformation, for alternative extensions of spin networks see |46j . where 
'ternary trees' are introduced to represent graphically the basic features of 'elliptic' 
coordinate sets on S 2 and S 3 (S denotes the standard d-sphere embedded in 
the (d + l)-dimensional Euclidean space) and of the corresponding harmonics. 
Interestingly, moving continuously along the edges of the abacus of Fig. Q] can be 
associated to the variation of the modulus of elliptic functions [33 133 132] • 

In Sec. 14.21 we have outlined how the concepts of reference frame transforma- 
tions and of alternative angular momentum coupling schemes in quantum mechan- 

automaton in computer science is a graph whose nodes encode 'internal states' 
while a link between two nodes represents an admissibile 'transition' between the corre- 
sponding states. 

2 A topological invariant is a quantity -typically a number or a polynomial- that de- 
pends only on the global topology of the geometrical object and not on its local metric 
properties. 



10 



ics lead to different representations of the quantum scattering matrix and provide 
a powerful guide for the analysis of atomic and molecular collisions. In particu- 
lar, we have exemplified atomic and molecular elastic and inelastic collisions, but 
extensions to reactive scattering are most interesting and extensive applications 
have been worked out. Dynamical calculations for the system He+H^ |50[ I51| and 
for the benchmark reaction F+H2 |52[ 153] have been performed, also including 
fine-structure and isotopic effects on reactivity. 

Also, a new class of entrance and exit channel indices in the scattering matrix 
has been worked out. Through the hyper spherical coordinate formulation referred 
to in Sec. I4.2[ the hyperrradial problem is essentially equivalent to that of scatter- 
ing from anisotropic potentials, and such a "stereodirected representation" of the 
scattering matrix can be used to derive information about the stereodynamincs 
of an atom-diatom reaction. A quantity that can be reconducted to properties 
measured in beam experiments on oriented molecules is the reaction probabilitiy 
as a function of the steric quantum number (see e.g. |26j). 

Since present quantum-mechanical calculations are becoming feasible for reac- 
tive encounters on realistic potential energy surfaces (the hyperquatization algo- 
rithm provides an efficient machinery in this direction) , stereodynamical properties 
exploiting the stereodirected representation have been reported, specifically for the 
reaction of HF with Li [54, 55J and of F with H2 |56j. 

Further perspectives concerning the extension of angular momentum theory to 
hyperspaces and the use of modern advances in the theory of orthogonal polynomi- 
als of a discrete variable have been reviewed |57j . Among applications, it is worth 
mentioning the possibility of representing polarization parameters by 'discrete' 
multipole moments [58l [59] and potential energy surfaces by orthogonal discrete 
basis sets [60] , 

The study of asymptotic expansions of 3nj symbols for n > 2, as well as of 
different types of asymptotics (in which only a few variables are large, while the 
other ones are kept 'quantized') represents a major challenge not only in the frame- 
work of the formal theory of hyper geometric polynomials and related hierarchies 
[I], but also in view of applications to specific physical problems arising in connec- 
tions with all the topics discussed in the previous Sections H] and A key example 
is the crucial occurrence of 9j symbols in the many center problem in quantum 
chemistry, either in Sturmian orbital or in momentum space approaches [61, .62] . 

Finally, for what concerns the issues discussed in Sec. [5l further algorithmic 
problems regarding (spin network-type) quantum geometry, topological quantum 
field theories in dimension 3 and associated 2-dimensional lattice models (as well 
as relations among them) are currently addressed [63] . 
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